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Abstract 
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We report a measurement of the time-dependent ratio of D° — > K + ir~ to 
Z)° — > K~ir + decay rates in D* + -tagged events using l.Ofb -1 of integrated lu- 
minosity recorded by the LHCb experiment. We measure the mixing parameters 
^ x 12 = (-0.9 ± 1.3) x 10" 4 , y' = (7.2 ± 2.4) x 10~ 3 and the ratio of doubly-Cabibbo- 

suppressed to Cabibbo-favored decay rates Rd = (3.52 ± 0.15) x 10~ 3 , where the 
uncertainties include statistical and systematic sources. The result excludes the 
no- mixing hypothesis with a probability corresponding to 9.1 standard deviations and 
represents the first observation of D° — D° oscillations from a single measurement. 
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Meson-antimeson oscillations are a manifestation of flavor changing neutral currents 
that occur because the flavor eigenstates differ from the physical mass eigenstates of 
the meson-antimeson system. Short-range quark-level transitions as well as long-range 
processes contribute to this phenomenon. The former are governed by loops in which 
virtual heavy particles are exchanged making the study of flavor oscillations an attractive 
area to search for physics beyond the standard model (SM). Oscillations have been observed 
in the K° — K° jl|, B° — B° [2| and B Q S — B° s [3) systems, all with rates in agreement with 
SM expectations. Evidence of D° — D° oscillations has been reported by three experiments 
using different D° decay channels [4]-[7]. Only the combination of these measurements 
provides confirmation of D° — D° oscillations, also referred to as charm mixing, with more 
than 5cr significance |8). While it is accepted that charm mixing occurs, a clear observation 
of the phenomenon from a single measurement is needed to establish it conclusively. 

Charm mixing is characterized by two parameters: the mass and decay width differences, 
Am and Ar, between the two mass eigenstates expressed in terms of the dimensionless 
quantities x = Am/T and y = Ar/2r, where T is the average D° decay width. The charm 
mixing rate is expected to be small, with predicted values of \x\, \y\ < 0(1CT 2 ), including 
significant contributions from non-perturbative long-range processes that compete with 
the short-range electroweak loops [9-12 . This makes the mixing parameters difficult to 



calculate and complicates the unambiguous identification of potential non-SM contributions 



in the experimental measurements 13 -15 . Current measurements give values towards the 
upper end of the range of theoretical predictions (8,13 . Together with the recent evidence 
of CP violation in charm decays 16, 17 , this can be interpreted as a hint for the presence 
of non-SM physics. 

In the analysis described in this Letter, D° — D° oscillations are observed by studying 
the time-dependent ratio of D° — > K + ir~ to D° — > K~ir + decay ratesQ The D° flavor at 
production time is determined using the charge of the soft (low-momentum) pion, 7r+, in 
the strong D* + — > D°-k+ decay. The D* + — > _D°(— » i^~7r + )7r+ process is referred to as 
right-sign (RS), whereas the D* + — > D°(—> K + -k~)ti+ is designated as wrong-sign (WS). 
The RS process is dominated by a Cabibbo-favored (CF) decay amplitude, whereas the 
WS amplitude includes contributions from both the doubly-Cabibbo-suppressed (DCS) 
D° K + -k~ decay, as well as D° — D° mixing followed by the favored D° — > K + n~ 
decay. In the limit of small mixing (\x\, \y\ <C 1), and assuming negligible CP violation, 
the time-dependent ratio, R, of WS to RS decay rates is approximated by [9] 



R(t)*R D + ^y> t - + X ^^(^\ , (1) 

where t/r is the decay time expressed in units of the average D° lifetime r, Rp is the ratio 
of DCS to CF decay rates, x' = x cos 5 + y sin 6, y' = y cos 6 — x sin 6, and 6 is the strong 
phase difference between the DCS and CF amplitudes. 

The analysis is based on a data sample corresponding to l.Ofb -1 of = 7 TeV pp 



collisions recorded by LHCb during 2011. The LHCb detector 18 is a single-arm forward 



lr rhe inclusion of charge-conjugated modes is implied throughout this Letter. 
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spectrometer covering the pseudorapidity range 2 < r] < 5, designed for the study of 
particles containing b or c quarks. Detector components particularly relevant for this 
analysis are the silicon Vertex Locator (VELO), which provides identification of displaced, 
secondary vertices of b- and c-hadron decays; the tracking system, which measures charged 
particles with momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% at 
100 GeV/c, corresponding to a typical mass resolution of approximately 8 MeV/c 2 for a 
two-body charm-meson decay; and the Ring Imaging Cherenkov (RICH) detectors, which 
provide kaon-pion discrimination. 

Events are triggered by signatures consistent with a hadronic charm decay. The 
hardware trigger demands a hadronic energy deposition with a transverse component of at 
least 3 GeV. Subsequent software-based triggers require two oppositely-charged tracks to 
form a D° candidate with a decay vertex well separated from the associated primary pp 
collision vertex (PV). Additional requirements on the quality of the online-reconstructed 
tracks, their transverse momenta (pr) and their impact parameters (IP), defined as the 
distance of closest approach of the reconstructed trajectory to the PV, are applied in the 
final stage of the software trigger. For the offline analysis only D° candidates selected by 
this trigger algorithm are considered. 

Selected D° candidates are combined with a track with p T > 300 MeV/c and p > 
1.5 GeV/c to form a D* + candidate. The D° daughter particles are both required to have 
p T > 800 MeV/c, p > 5 GeV/c and x 2 (IP) > 9. The x 2 (IP) is defined as the difference 
between the \ 2 of the PV reconstructed with and without the considered particle, and 
is a measure of consistency with the hypothesis that the particle originates from the 
PV. The D° candidate is also required to have pr > 3.5 GeV/c. Contamination from D 
mesons originating from 6-hadron decays (secondary D) is reduced by requiring the x 2 (IP) 
of the D° and of 7r+ candidates to be smaller than 9 and 25, respectively. In addition, 
the RICH system is used to distinguish between pions and kaons and to suppress the 
contamination from misidentified two-body charm decays in the sample. Backgrounds 
from misidentified singly Cabibbo-suppressed decays are specifically removed by requiring 
the D° candidate mass reconstructed under the K + K~ and 7r + 7r~ hypotheses to differ by 
more than 40 MeV/c 2 from the known D° mass [l9j. Contamination from electrons to the 
soft pion sample is also suppressed using particle identification information. Finally, it is 
required that the D° and the 7r+ form a vertex, which is constrained to the measured PV. 
Only candidates with reconstructed Ktt mass within 24 MeV/c 2 of the known D° mass 
and with reconstructed D°tt^ mass below 2.02 GeV/c 2 are considered further. The D°ir^ 
mass, M(_D°7r+), is calculated using the vector sum of the momenta of the three charged 
particles and the known D° and ir + masses [l9]; no mass hypotheses for the D° daughters 
enter the calculation, ensuring that all two-body signal decays have the same M(Z)°7r+) 
distribution |20|. Events with multiple RS or WS D* + candidates occur about 2.5% of the 
time, and all candidates are kept. 

Figure [l] shows the M(D°n+) distribution for the selected RS and WS candidates. 
Overlaid is the result of a binned x 2 fit used to separate the D* + signal component, with a 
mass resolution of about 0.3 MeV/c 2 , from the background component, which is dominated 
by associations of real D° decays and random pions. The signal mass shape is modelled as 
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Figure 1: Time- integrated -D°vr+ mass distributions for the selected RS D° — » fC~7r + (left) 
and WS L> — > K + tt~ (right) candidates with fit projections overlaid. The bottom plots 
show the normalized residuals between the data points and the fits. 



the sum of one Johnson Sy [21] and three Gaussian distributions, which account for the 
asymmetric tails and the central core of the distribution, respectively. The background is 
described by an empirical function of the form [M(D°tt+) — m ] a e ~ b [ M ^ D ° ns )~ m o]^ w h ere 



the threshold mo is fixed to the sum of the known D° and tt + masses 19 . We reconstruct 
approximately 3.6 x 10 4 WS and 8.4 x 10 6 RS decays. To determine the time-dependent 
WS/RS ratio the data are divided into thirteen D° decay time bins, chosen to have a 
similar number of candidates in each bin. The decay time is estimated from the distance L 
between the PV and the D° decay vertex and from the D° momentum as t/r = m D oL/pr, 



where m^o and r are the known D° mass and lifetime 19 , respectively. The signal 
yields for the RS and WS samples are determined in each decay time bin using fits to 
the M(D°7r+) distribution. The shape parameters and the yields of the two components, 
signal and random pion background, are left free to vary in the different decay time bins. 
We further assume that the M(D°n^) signal shape for RS and WS decay are the same 
and therefore perform first a fit to the abundant and cleaner RS sample to determine the 
signal shape and yield, and then use those shape parameters with fixed values when fitting 
for the WS signal yield. The signal yields from the thirteen bins are used to calculate the 
WS/RS ratios, shown in Fig. [2J and the mixing parameters are determined in a binned x 2 
fit to the time-dependence according to Eq. Q. 

Since WS and RS events are expected to have the same decay-time acceptance and 
M(D°7r+) distributions, most systematic uncertainties affecting the determination of 
the signal yields as a function of decay time cancel in the ratio between WS and RS 
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events. Residual biases from non-cancelling instrumental and production effects, such as 
asymmetries in detection efficiencies or in production, are found to modify the WS/RS 
ratio only by a relative fraction of O(10~ 4 ) and are neglected. Uncertainties in the distance 
between VELO sensors can lead to a bias of the decay-time scale. The effect has been 



estimated to be less than 0.1% of the measured time 22 and translates into relative 
systematic biases of 0.1% and 0.2% on y' and x' 2 , respectively. At the current level of 
statistical precision, such small effects are negligible. 

The main sources of systematic uncertainty are those which could alter the observed 
decay-time dependence of the WS/RS ratio. Two such sources have been identified: (1) 
secondary D mesons, and (2) backgrounds from charm decays reconstructed with the 
wrong particle identification assignments, which peak in M(D°7r^) and are not accounted 
for in our mass fit. These effects, discussed below, are expected to depend on the true 
value of the mixing parameters and are accounted for in the time-dependent fit. 

The contamination of charm mesons produced in 6-hadron decays could bias the time- 
dependent measurement, as the reconstructed decay time is calculated with respect to 
the PV, which, in this case, does not coincide with the D° production vertex. When this 
secondary component is not subtracted, the measured WS/RS ratio can be written as 
R(t) [1 — A B (t)}, where R(t) is the ratio of promptly- produced candidates according to 
Eq. ([!]), and A B (t) is a time-dependent bias due to the secondary contamination. Since 
R(t) is measured to be monotonically non-decreasing [8] and the reconstructed decay time 
for secondary decays overestimates the true decay time of the D° meson, it is possible to 
bound A B (t), for all decay times, as 



< A B (t) ^ /| s (t) 



R(t)\ 



(2) 



where f B s (t) is the fraction of secondary decays in the RS sample at decay time t. The lower 
bound in Eq. ^ corresponds to the case when the parent 6-hadron decays instantaneously 
and the reconstructed D° decay time is the true decay time. The upper bound corresponds 
to the case when the D° decays instantaneously and the reconstructed decay time t is 
entirely due to the 6-hadron lifetime. Since A B ^ 0, it follows that the background from 
secondary D decays decreases the observable mixing effect. To include the corresponding 
systematic uncertainty, we modify the fitting function for the mixing hypothesis assuming 
the largest possible bias from Eq. ^ . A measurement of the secondary fraction is done by 
fitting the x 2 (IP) distribution of the RS D° candidates in bins of decay time. The value of 
f B s (t) is constrained in the time-dependent fit to the measured fraction, which increases 
almost linearly with decay time from (0.0 ± 0.5)% up to (14 ± 5)%, for a time-integrated 
value of (2.7 ±0.2)%. 

Background from incorrectly reconstructed D meson decays, peaking in the M(D°n^) 
distribution, arises from D* + decays for which the correct soft pion is found but the D° is 
partially reconstructed or misidentified. This background is suppressed by the use of tight 
particle identification and two-body mass requirements. From studies of the events in the 
D° mass sidebands, we find that the dominant peaking background is from RS events that 
survive the requirements of the WS selection; they are estimated to constitute (0.4 ± 0.2)% 
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Figure 2: Decay-time evolution of the ratio, R, of WS D° — >■ if + 7r~ to RS -D — >■ K~n + 
yields (points) with the projection of the mixing allowed (solid line) and no-mixing (dashed 
line) fits overlaid. 

of the WS signal. This contamination is expected to be independent of decay time and, if 
neglected, would lead to a small increase in the measured value of Rp. From the events 
in the D° mass sidebands, we derive a bound on the possible time dependence of this 
background, which is included in the fit in a similar manner to the secondary background. 
The x 2 t na t is minimized in the fit to the WS/RS decay-time dependence is 

X (ri,ti,(Ti\0) = I J + X B \ e ) +Xp{0), (3) 

where and crj are the measured WS/RS ratio and its statistical uncertainty in the decay 
time bin i, respectively. The decay time ti is the average value in each bin of the RS 
sample. The fit parameters, 0, include the three mixing parameters (Rd, y', x' 2 ) and five 
nuisance parameters used to describe the decay time evolution of the secondary D fraction 
(As) and of the peaking background (A p ). The nuisance parameters are constrained to the 
measured values by the additional Xb an d Xp terms, which also includes their correlations. 

The analysis procedure is defined prior to fitting the data for the mixing parameters. 
Measurements on pseudo-experiments that mimic the experimental conditions of the data, 
and where D° — D° oscillations are simulated, indicate that the fit procedure is stable and 
free of any bias. 

The fit to the decay-time evolution of the WS/RS ratio is shown in Fig. [2] (solid 
line), with the values and uncertainties of the parameters Rd, y' and x' 2 listed in Table [l] 
As the dominant systematic uncertainties are treated within the fit procedure (all other 
systematic effects are negligible), the quoted errors account for systematic as well as 
statistical uncertainties. When the systematic biases are not included in the fit, the 
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Table 1: Results of the time- dependent fit to the data. The uncertainties include statistical 
and systematic sources; ndf indicates the number of degrees of freedom. 



Fit type 


Parameter 


Fit result 


Correlation coefficient 


( X 2 /ndf) 




(io- 3 ) 


Rd y' x' 2 


Mixing 


Rd 


3.52 ±0.15 


1 -0.954 ±0.882 


(9.5/10) 


y' 


7.2 ±2.4 


1 -0.973 




x' 2 


-0.09 ±0.13 


1 


No mixing 


Rd 


4.25 ± 0.04 




(98.1/12) 
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Figure 3: Estimated confidence-level (CL) regions in the (x' 2 ,y') plane for 1 — CL = 0.317 
(lcr), 2.7 x 10~ 3 (3<r) and 5.73 x 10 -7 (5c). Systematic uncertainties are included. The 
cross indicates the no- mixing point. 



estimated uncertainties on R D , y' and x' 2 become respectively 6%, 10% and 11% smaller, 
showing that the quoted uncertainties are dominated by their statistical component. To 
evaluate the significance of this mixing result we determine the change in the fit x 2 when 
the data are described under the assumption of the no-mixing hypothesis (dashed line 
in Fig. [2]). Under the assumption that the x 2 difference, A% 2 , follows a x 2 distribution 
for two degrees of freedom, Ax 2 = 88.6 corresponds to a p-value of 5.7 x 10~ 20 , which 
excludes the no-mixing hypothesis at 9.1 standard deviations. This is illustrated in Fig. [3] 
where the lcr, 3cr and 5a confidence regions for x' 2 and y' are shown. 

As additional cross-checks, we perform the measurement in statistically independent 
sub-samples of the data, selected according to different data-taking periods, and find 
compatible results. We also use alternative decay-time binning schemes or alternative 
fit methods to separate signal and background, and find no significant variations in the 
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estimated parameters. 

In conclusion, we measure the decay time dependence of the ratio between D° — » K + n~ 
and D° — > K~tt + decays using 1.0 fb~ of data and exclude the no-mixing hypothesis at 
9. 1 standard deviations. This is the first observation of D° — D° oscillations in a single 
measurement. The measured values of the mixing parameters are compatible with and 



have substantially better precision than those from previous measurements 4 6 , 23 
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